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A  CLASS  OF  IWO -PHASE  STATE  EQUATIONS 


ABSTRACT 

The  properties  of  a  particular  class  of  elementary,  two-phase 
state  equations  are  found  to  be  simply  related.  The  class  includes 
some  familiar  equations,  such  as  the  Van  der  Waals,  Berthelot,  and 
Clausius  types,  as  well  as  a  form  recently  suggested  by  J.  J.  Martin. 
Conditions  for  stable  equilibrium  are  given,  and  some  examples  are 
compared.  A  table  of  values  is  included  which  can  be  used  with 
elementary  transformations  to  find  the  equilibrium  coexistence 
properties  of  any  member  of  the  set.  A  drawback  for  this  type 
equation  along  the  critical  isotherm  is  noted. 
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I.  INTRODUCTION 


There  exists  a  substantial  number  of  elementary  pressure  state 
equations  which  predict  a  liquid-vapor  phase  transition  and  a  critical 
point.  It  is  the  purpose  of  this  paper  to  show  that  many  of  these, 
though  not  all,  belong  to  a  single  class,  and  that  the  various  thermo¬ 
dynamic  predictions  between  class  members  are  simply  related.  Some 
familiar  equations  belong  to  this  group,  including  the  Van  der 
Waal s -Maxwell,  [V-M],  equation,  which  will  serve  as  the  basis  for 
comparison.  A  direct  method  can  be  demonstrated  relating  the  limiting 
values  of  thermodynamic  quantities  in  the  critical  region. 

The  form  of  equation  treated  here  is 

P(V,T)  =  (RT/[V-b(T)])  -  a(T)/[V+c(T)]2  (l) 

where  V  ^  b  (T)  >  -  c(t),  and  T  >  0;  P,  V,  and  T  are  pressure,  specific 
volume  and  temperature,  respectively.  R  is  the  gas  constant,  and  a(T), 
b(T)  and  c(t)  are  analytic  functions  of  positive  T.  Some  restrictions 
are  needed  on  the  latter  functions.  For  instance,  b(T)  is  assumed  to 
be  positive  and  larger  than  -c(T);  and  a(T)  must  be  positive  for  all  T 
less  than  some  value  exceeding  the  critical  temperature,  Tc ,  in  order 
to  produce  a  critical  point.  It  is  shown  in  the  Appendix  that  b(T) 
must  be  a  constant,  which  is  henceforth  assumed. 

Eq.  (1)  is  expressed  in  reduced  form  by  use  of  the  scaled 
variables  3  =  P/Pc ,  y  =  T/Tc  ,  T]  =  l/a  =  V/Vc  (this  notation  corresponds 
to  that  of  Reference  [1],  with  subscript  denoting  critical  value)  as 

eoi,Y)  =  (i/Zc)  f (y/m-Bj)  -  (9/8)A(Y)/m+c(Y)]2}  (2) 

where  T|  s  B  >  -  C(y)  and  y  >  0,  and  Zc  =  PCVC/(RTC)  is  the  critical  com¬ 
pressibility.  The  transformations  from  a,  b,  and  c  to  A,  B,  and  C  are 
(9/8)  A(y)  =a(T)/[RTcVc ],  B=b/Vc  and  C(y)=c(T)/Vc .  It  follows  that 
A(y)  is  positive  for  a  range  of  y  larger  than  unity.  The  variables 
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P,  Y,  T|  and  a  all  are  unity  at  the  critical  point.  It  is  seen  that 
Eq.  (l)  approximates  the  ideal  gas  equation  when  V  becomes  very  large. 


In  order  that  the  state  of  the  system  be  stable,  certain  relations 
must  hold  for  the  thermodynamic  functions2* .  These  are:  Cv  >  0,  where 
Cv  is  the  specific  heat  at  constant  volume;  and  (,=7^,  <  0,  except  at 
the  critical  point.  At  the  critical  point,  the  conditions  become 


(i)v  -  0  - 


ft*  ft 


Imposing  these  latter  conditions  at  the  critical  point,  one 
obtains  the  relations: 


Ai  =  1+Cij  Bi  =  (l/3)(l-2Ci);  and  Zc  =  (3/8)/(l+Cx  )  (3) 

where  a  subscript  i_  on  A,  B,  or  C  denotes  the  value  at  the  critical 
temperature,  y  =  1.  Thus,  the  function  C(y)  fixes  the  critical 
compressibility,  Zc . 

9  B 

One  observes  that  for  y  <  1,  there  is  a  domain  of  (T| ,Y )  where  gyj- 
from  Eq.  (2)  is  positive,  contrary  to  the  stability  requirement. 

Thus,  Eq.  (2)  cannot  represent  a  stable  thermodynamic  pressure  function 
in  this  domain;  instead,  the  vapor  pressure  function,  P*(y),  independ¬ 
ent  of  T|,  replaces  Eq.  (2)  as  the  pressure  function  here.  The  usual 
method  of  defining  PA  (y)  so  it  is  continuous  with  the  P(T|,y)  of  Eq.  (2) 
is  by  means  of  the  Maxwell  rule:  for  any  fixed  y  £  1, 

My) 

3Cn,y)cfn  =  Pa(y)CMy)  -  My)},  (4) 

My) 

with 


e(My)  =  0(My)  =  My)  • 


^References  are  listed  on  page  27 . 
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The  quantities  Tl3  (y)  and  %  (y)  are  points  on  the  saturated  liquid  and 
vapor  loci,  respectively.  For  any  fixed  y  ^  1,  the  domain  ^3  ^  ^  % 

is  the  region  where  liquid  and  vapor  coexist  at  the  same  pressure, 
temperature,  and  chemical  potential.  The  boundaries  of  the  coexistence 
region  are  also  defined  by  Eq.  (4).  The  use  of  the  Maxwell  rule 
implies  that  the  Helmholtz  free  energy  is  a  unique  analytic  function  of 
V  and  T  where  Eq.  (l)  is  the  proper  pressure  function3.  With  this 
analyticity,  Eq.  (4)  insures  that  the  pressure  and  chemical  potential 
are  continuous  functions  of  V  and  T  across  the  coexistence  boundaries. 
To  avoid  confusion,  a  subscript  *_  will  hereafter  be  used  on  a  thermo¬ 
dynamic  variable  within  the  coexistence  region,  while  a  subscript  b_ 
will  refer  to  the  value  outside  coexistence. 

II.  TRANSFORMATIONS 


For  compactness,  let 

f(v)  =  3[b-k:(y)  ]/A(y)  .  (5) 

For  y  >  1,  must  never  be  positive,  or  the  stability  requirement  is 
violated.  From  Eq.  (2), 

Ze  ^  -  (V/[H-B]>){1-G(H,Y)},  (6) 

where 

gP1,y)  =  (9/1*)  A(Y)[H-Bf/(Y[T1->C]3)  .  (7) 

Since  both  Zc  and  y  are  positive,  Eq.  (6)  shows  is  negative, 
providing  G(T|,y)  is  less  than  1  for  any  T]  £  B.  This  latter  condition 
is  easily  satisfied  when  A(y)  is  negative;  but  A(y)  was  previously 
required  to  be  positive  for  all  y  less  than  some  value  y0  exceeding 
one.  In  this  range  of  y,  i.e.,  (0  <  y  <  y0),  G(T1,y)  is  positive  and  has 
a  maximum  value  with  respect  to  T],  GM(y),  at  Tj  =  3B+2C(y),  where 
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Gm(y)  =  A(y)/[3y(B+C)]  =  l/[Yf(Y)]. 


(8) 


Forcing  Yf  (y)  >  1  for  Y0  >  Y  >  1  will  insure  that  is  non-positive 

For  y  <  lj  one  wishes  to  satisfy  Eq.  (4)  in  the  CHjY)  domain  of 
coexistence.  The  mathematical  function  j'Y)  of  Eq.  (2)  must  there 
fore  have  zeros  in  this  domain,  so  Yf(Y)  <  1  for  y  <  1.  At  Y  =  1, 

Yf  (y)  is  unity,  as  may  be  seen  from  the  relations  already  given  for 
Y  =  1. 

Defining 


x(y)  s  (3/8)/ [Zc  f2 (y)a(y)  ] , 

(9) 

and  the  transformations 

u(Y)  s  Yf  (y)  5 

(10) 

v(ii,y)  5  rn+c (y ) 3/ [f (y )a(y ) ]> 

(11) 

w(t1,y)  s  3(t1jy)/x(y) j 

(12) 

wa(y)  =3a(y)/x(y). 

(13) 

with  the  convention  that  w(i|,Y)  has  a  subscript  b_  when  3(^1, y)  does, 
one  finds  that  Eq.  (2)  transforms  into 

w(v,u)  =  f8u/(3v-l)}  -  3/v2  .  (l4) 

Moreover,  for  any  fixed  u  ^  1,  Eq.  (4)  transforms  into 

vx  (u) 

w(v,u)dv 

t 

v3(u) 

with  w(va  ,u)  =  Wg  (vj  ,u)  =  w(v3  ,u)  =  wB  (v3  ,u)  =  wA  (u). 


=  wA  (u)  fv^u)  -  v3(u)}  (15) 
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Eq.  (l4)  has  exactly  the  form  of  the  reduced  Van  der  Waals  pressure 
function,  providing  one  interprets  w,  v,  and  u  as  reduced  pressure, 
volume,  and  temperature  respectively.  All  the  functions  wB  ,  wA  ,  u,  and 
v,  are  unity  at  the  critical  point  where  "0  =  1  and  y  =  1. 

The  Jacobian,  J  =  3  transformations  of  the  independent 

variables  y  and  4  to  the  new  independent  variables  u  and  v  reduces  to 
because  u  is  independent  of  0 .  Using  Eqs.  (lO)  and  (ll)  this 
Jacobian  is 

j  =  Ei+v  g/f(v))/A(y),  (16) 

and  where  J  or  J"1  is  neither  zero  nor  infinite,  the  mapping  is  single 
valued.  It  will  be  assumed  in  the  further  development  that  for  y  ^  1 
this  mapping  is  single  valued,  so  the  necessary  conditions  on  A(y),  B 
and  C(y)  to  accomplish  this  purpose  will  be  assumed  henceforth.  The 
thermodynamic  properties  of  any  state  equation  of  this  class  for  y  £  1 
can  now  be  expressed  in  terms  of  those  corresponding  properties  known 
for  the  [V-M]  system,  and  the  given  transformation  Eqs.  (10-13). 

The  method  for  finding  the  properties  in  the  coexistence  region, 
including  their  limits  at  the  critical  point,  is  that  demonstrated  by 
Barieau" .  For  fixed  u  £  1,  Eq.  (l4)  becomes  a  cubic  polynomial  in  v; 
one  finds  from  it  the  largest  and  smallest  values  of  v  (vi (u)  and 
v3 (u) ,  respectively)  which  simultaneously  satisfy  Eq.  (15).  Letting 
|  =  l/v,  one  then  obtains,  for  u  ^  1, 

wa(5i,§3)  =  5i5s(3-?s-5i)»  (17) 

8u(?i,?3)  =  (li+?3 )(3-^i )(3-?3 ).  (18) 

Differentiating  both  sides  of  Eq.  (15)  by  u  results  in  the  differential 
equation: 


U  ^  =  WA  +  35x5a  =  5i5a(6-e3-?x); 


(19) 
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similarly  differentiating  Eqs.  (17)  and  (l8),  and  using  Eq.  (17)  gives 
one 


(20) 

u  It3-  =  ?a(3-5a)/(3-2?3-?i).  (21) 


The  relation 


u 


d2w* 


du 


du  ’ 


(22) 


follows  directly  from  Eq.  (19).  Eqs.  (17-22)  have  the  same  form  as 
the  corresponding  ones  of  Reference  [1]  used  to  display  properties  of 
the  [V-M]  fluid,  providing  wA ,  u,  and  §  here  are  replaced  by  0,  Y,  and 
or  to  conform  with  the  notation  there . 

Proceeding  analogously,  and  defining  y  =  -  1  and  x  =  1  -  ?i  , 

one  can  use  Eqs.  (l4),  (15),  (17)  and  (18)  to  obtain  y  as  a  function  of 
x;  this  function  can  be  expanded  in  the  power  series 

00 

y=^  aiX1  ,  (23) 

i=l 


where  Barieau1  finds  ai=l,  a^- 1/5?  and  a3 -1/ 25 .  The  value  of  a4  was 
subsequently  evaluated  by  the  present  author  to  be  19/350. 

All  quantities  5i  ,  ?3  ,  wA  ,  etc.,  (as  functions  of  u)  have  been 
tabulated  previously  for  the  [V-M]  system.  They  may  be  used  to  find 
limiting  thermodynamic  values  in  the  critical  region  for  any  equation 
of  the  class  considered  here;  a  brief  table  is  in  the  Appendix. 

Some  inverse  transformations  are 

\  =  1/oTj  =  Vj  f  (y)  A(y )  -  C(y)  ,  (j=l,3)  (24) 
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PbC^jY)  =w9(v,u)x(y) 


(25) 


(y)  =  wA (u)X(y) 


§ 


h  (t1?y 


W 


(26) 

(27) 


■dPi  (V.) 
dY 


dX 

dY  • 


(28) 


From  Eqs .  (27)  and  (28),  it  can  be  seen  that  {gy6  is  e(lual  to 

>  because  it  is  known  from  the  Van  der  Waals  function  that  at  the 
cry 

point  (v,u)  =  (l,l):  (a),  wB  and  wA  are  each  equal  to  one;  (b), 

and  '(d)  are  each  equal  to  four;  and  (c),  is  zero. 

Therefore,  for  any  equation  of  this  class,  the  temperature  derivative 
of  the  vapor  pressure  function  is  continuous  at  the  critical  point  with 
that  of  the  critical  isometric,  in  reduced  variables.  This  continuity 
does  not  necessarily  hold  for  higher  order  derivatives. 


Ill  EXAMPLES  AND  COMPARISONS 

When  A(y),  B  and  C(y)  are  specified,  one  can  recognize  some 
well-known  state  equations.  Some  examples  are  shown  in  Table  I. 

Example  6,  with  A(y)  a  linear  function  of  y,  comes  from  an  equation 
recently  suggested  by  Martin4*.  When  C(y)  is  a  constant,  as  is  true 
for  all  equations  in  Table  I,  one  can  use  the  definitions  to  show  that 
the  product  f(y)x(y)  is  unity.  In  the  following,  certain  thermodynamic 
predictions  for  some  of  these  examples  will  be  compared. 

A.  Coexistence  Boundary. 

A  comparison  is  made  in  Figure  1  of  the  reduced  density  on  the 
coexistence  boundary  for  the  state  equations  numbered  2,  4,  and  6  in 


15 


Table  I.  Example  state  equations  and  transformations 


a.  Shown  in  the  form  P(V,T)  =  [RT/(V-b)]  -  a(T)/(V+c)2,  with 
constants  aj  ,  bj  and  Cj  for  the  number  (j)  equation.  For  all  shown 
(except  No.  (l),  where  both  b  and  a(T)  are  zero},  bj  =  Vc  [l-l/(4zc)] 
IX-  and  a  =  Vc[-l+3/(8Zc)];  n,  =  (9/64)(RTc  )2/Pc  • 


No. 

(J) 

Name 

p(v,t) 

aJ 

zc 

(1) 

Ideal  Gas 

rt/v 

-  Unde 

fined- 

(2) 

Van  der  Waals 

[RT/(V-tfe  )  ]-ae  /V2 

3  ng 

*  3/8 

(3) 

(translated) 
Van  der  Waal  s 

[RT/(V-^  )]-a3/(V+c3)2 

3  n3 

>  1/4 

(4) 

Berthe lot 

[RT/ (V-b4  )  ]-&*/ (TV2  ) 

3  r^Tc 

=  3/8 

(5) 

Clausius 

[RT/(V-bB  )  ]-%  /  [T(V+cb  )2 

3  n^Tc 

>  1/4 

(6) 

Mart in4  (A) 

(RT/ (V-be )  ]-ae (4  Tc -T)/(V+c6 )2 

ne/Tc 

>  1/4 

(7) 

Martin  (b) 

[RT/(V-b7)]-a7(5  Tc-2T)/(V+c7)2 

n7/Tc 

>  1/4 

b.  Same  state  equations  (omitting  No.  (l)}  in  the  reduced  form 
P  =  (l/Zc)([Y/(1fl-B)]-(9/8)A(Y^'n^c)2},  where  P  =  P/Pc  , 

Y  =  T/Tc  and  T]  =  T/Tc  .  For  all  shown  here,  B  =  1  -  l/(4Zc), 
C  =  -  1  +  3/(8zc)  (independent  of  y! ,  and  X(y)  =  l/f(Y). 


No. 

(j) 

zca(y) 

Zc 

f(Y) 

(2b)  Lftb)\ 

'  dY  /c  !  V 

(2) 

3/8 

=  3/8 

1 

4 

(3) 

3/8 

>  1/4 

1 

4 

(4) 

3/(8y) 

=  3/8 

Y 

7 

(5) 

3/(8y) 

>  1/4 

Y 

7 

(6) 

(4-y)/8 

>  1/4 

3/(4-y) 

5 

(7) 

(5-2y)/8 

>  1/4 

3/(5-2y) 

6 

16 


HOI 
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Figure  1.  Reduced  temperature  vs.  saturation  densities.  Numbers  refer  to  state 
equations  in  Table  I.  For  (6),  the  value  of  Zc  is  assumed  to  be  l/3. 


Table  I,  using  the  given  functions  with  Eq.  (24);  similar  plots  for 
examples  3  and  5  of  Table  I  can  easily  be  obtained  from  those  of  2  and 
4  respectively,  when  the  value  of  Ci  is  assigned. 

B.  Vapor  Pressure,  (y  ^  l). 

Eq.  (26)  gives  the  vapor  pressure,  PA(y),  as  a  function  of  wA(u), 
the  [V-M]  vapor  pressure  at  the  reduced  temperature  u,  and  the  specific 
transforming  function  X(y).  These  are  shown  in  Figure  2.  For  state 
equations  having  the  same  f(y)  functions  and  the  same  X(y)  functions, 
the  vapor  pressure  is  the  same  function  of  y.  Thus,  the  state 
equations  numbered  2  and  3  in  Table  I  have  the  same  vapor  pressure,  as 
do  those  numbered  4  and  5. 

C .  Heat  of  Vaporization. 

For  T  ^  Tc ,  the  latent  heat  of  vaporization  Aqa (t)  is  found  by 
using  the  Clapeyron  relation,  Aqa  (t)  =  T  [Vi  (t)  -  V3(t)];  in 
reduced  variables  this  becomes 

Aqa(t)/[rtc]  =  zc  y  1^-  Pli  (y)  -  “My)]  •  (29) 

Using  the  inverse  transformations, 

aq*  (t)/[rtc  ]  =  f aqa  (t)/  [rtc  ]}0  {(|~^)pj/fs  (y )} 

+  (uwA  (u)[vi  (u)  -  v3(u)]}{zcA  (30) 

where  the  first  factor  on  the  right  is  the  latent  heat  of  vaporization 
for  the  [V-M]  equation  at  the  reduced  temperature  u.  The  functions 
wA(u),  vi (u)  and  v3 (u)  are  also  known;  as  with  the  vapor  pressure 
function,  the  heat  of  vaporization  is  the  same  function  of  y  for  state 
equations  having  the  same  f(y)  and  X(y).  Comparisons  are  shown  in 
Figure  3;  the  values  for  the  Berthelot  and  Clausius  equations  become 
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r 

Figure  2.  Vapor  pressure  vs.  temperature  for  equations  numbered 
in  Table  I.  Vapor  pressure  is  independent  of  Zq  for 
the  examples  shown. 
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r 


Figure  3.  Latent  heat  of  vaporization  vs.  reduced  temperature 
for  equations  numbered  in  Table  I.  Latent  heat  is 
independent  of  Zc  for  examples  shown. 
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unrealistically  large  at  small  y.  Obviously,  Aqa (t)  vanishes  at  T  =  Tc 
for  all  equations.  The  entropy  of  vaporization,  ASa(t),  is  simply 
computed  from  Aqa  (t)  by  using  T  Asa(t)  =  Aqa(t). 


D.  Specific  Heat  at  Constant  Volume,  . 

With  E(V,T)  the  internal  energy  for  unit  mass  (here  taken  to  be 
the  molar  mass)  the  constant  volume  specific  heat  is  Cv  (v,T)  =  \^)  ? 

so  that  \^yL)  =  d'ydr  •  The  form  of  Cv  within  coexistence  is  Cva(V,T), 

and  it  differs  from  that  of  the  single  phase,  Cvb(V,T).  By 
differentiating  the  thermodynamic  relation 


(31) 


one  obtains 


53E 

dlBV 


=  T 


;  using  reduced  variables  and  equating  the 


above  mixed  derivatives,  one  gets: 


9T1 


(32) 


Performing  a  partial  integration  with  respect  to  T]?  one  obtains 


cv8(T],y)/r  =  Cy  (y)/r  +  YF(^,y) 


(33) 


where,  with  a  prime  hereafter  signifying 


F(T|,y)  =  f  (9/8)/[^+C(y)  ]3}{A"  [T)+C(y)  ]2  -[2A'C'+AC']pWc(Y)] 

+  2  A(C  /)3 }  .  (34) 


Here  Cv (t)  is  the  specific  heat  of  the  single  phase  at  very  low  density 
(T|  -»  oo,  f(T|,y)  -*  o]  ;  for  a  monatomic  vapor  at  moderate  temperature 
C°(t)  is  3R/2.  Since  CVB  must  always  be  positive  in  the  single  phase, 
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further  restrictions  are  imposed  on  the  permissible  functions  A(y)  and 


c(y). 

than 

from  Eq.  (34)  one  obtains, 


It  is  for  this  reason  that  B  is  forced  to  be  a  constant  rather 
than  a  function  of  y,  as  shown  in  the  Appendix.  Because  2^  =  (|t[)y 


Zc  =  f  (9/8)/D1+C(y)  ]4H-  A^H^Cy)]8  +  2t2A/C/tAC',]n^C(Y)] 

-6A(C')3}.  (35) 


■^2  Q 

If  C'(y)  and  A"(y)  are  identically  zero,  then  and  f(T1,y) 

are  identically  zero  in  the  single  phase. 

Within  coexistence,  T  ^  Tc ,  V3(t)  ^  Vi(t),  the  specific  heat  is 
Cva(V,T);  it  can  be  found  by  using  Ea(V,T)  and  Pa(V,T)  in  Eq.  (3l)5 
integrating  it  at  fixed  T  ^  Tc  from  V3 (t)  to  some  V  within  coexistence, 
and  differentiating  the  result  by  T.  Both  Ea(v,t)  and  CvA(V,T)  are 
linear  functions  of  V.  Using  the  fact  that  the  internal  energy  and 
pressure  functions  are  continuous  across  the  boundary  of  coexistence, 
one  can  obtain  in  reduced  variables 


(cva(v,t)  -  c5(t)}/r  =  1(T1,y),  (36) 

where 

rn,Y)=Y^3  '(y)  &a'(y)}  +  YF(Tl3,Y)+VPA//(Y)m-Tl3(Y)].  (37) 

In  general  is  not  equal  to  3  A  ' (y)  at  (^jY),  so  the  constant 

volume  specific  heat  exhibits  a  discontinuity  across  the  coexistence 
boundary. 

Recalling  that  the  are  functions  of  u,  one  can  use  the  inverse 
transformations  and  arrange  the  terms  to  obtain: 
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rOl>Y)  =  {r(v,u)}ok1  (y)  +  (3-§3)(?3-?1)fk8(Y)-(3-2§3-?1)k3(Y)}  (38) 


+  f  (?3-?)/(?3?)}fu  —A-  k*  (y)  +UWA  kg  (y)}  +yF2(?3>y). 

Here  the  ki  are  functions  of  Y,  and: 

\  =  fu'/f)3 ;  s  6zcu'c7x/f 

k3  s  SZcXfC'f/f^A];  k*  s  AZC  [Xu"  +  2X/u/] 

kg  =  AZCX";  and  with  m  =  §3/[fA], 

F8(53,y)  =  (9/8)  m  [A"- [2A'C'+AC/0  m  +  2A[C/]3m33. 

{r(v,u)}0  is  the  corresponding  quantity  for  the  [V-M]  system  at  the 
point  (v,u)  inside  coexistence.  The  value  of  r (T) , y)  at  the  critical 
point  is  easily  found  because  both  f  and  ^  become  equal  to  f3  there, 
and  the  second  and  third  terms  vanish.  The  critical  value  of 
fr(i,i)}0  is  known  to  be  4.5,  and  neither  k^ (y)  nor  F2 (l,l)  is  infinite 
there.  Therefore,  r (T| , y )  is  always  finite  at  the  critical  point,  and 
state  equations  of  this  class  cannot  show  an  infinite  CVA  here,  as 
suggested  by  some  experimental  data.  Along  the  saturated  liquid  locus 
(C  =?3 ) ,  the  third  term  of  Eq.  (38)  vanishes.  When  C7  =  0  (as  happens 
for  all  example  equations  in  Table  I),  Eq.  (38)  simplifies  along  the 
saturated  liquid  locus,  (T^y)?  to: 


I’Pls.Y)  =  0'(v3,u)}0{u7f}3+(9/8)  5a  (u)  {YA'7[fA]}.  (39) 


k.  Vaporizing  Signal  Speed 

The  speed  of  a  small  amplitude  adiabatic  expansion  wave  (with 
partial  vaporization)  propagating  into  the  saturated  liquid5  is  denoted 

TTJ  T_ 

by  aA  (V3  , T ) ,  and  is  equal  to  V3  {t/Cv  a  (V3  ,T)}2  ;  in  reduced  notation, 
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dPj 

dY 


(40) 


aA(V-3,T)/[RTc]2  =  ZcT13 


(y/rm 


3>y ) 


Cv/R]}1 


A  comparison  of  these  speeds  is  shown  in  Figure  4  for  the  sample 
equations  numbered  2,  4,  and  6  in  Table  I,  with  the  assumption  that 
C° /R  is  3/2.  The  limiting  values  at  the  critical  point  for  certain 
coexistence  properties  of  sample  state  equations  are  shown  in  Table  IX. 

IV.  DISCUSSION 

The  predicted  thermodynamic  properties  of  coexistence  for  this 
class  of  elementary  state  equations  are  seen  to  be  conveniently  done  in 
terms  of  the  known  [V-M]  system.  The  use  of  different  functions  A(y), 

B  and  C(y),  might  improve  the  agreement  between  prediction  and  obser¬ 
vation.  However,  this  class  of  state  equation  contains  one  serious 
drawback  in  the  critical  region  which  cannot  be  overcome  by  modifying 
these  functions.  This  defect  is  that  this  entire  class  predicts  |p-Pc| 
to  vary  as  the  third  power  of  |v-Vc|  on  the  critical  isotherm;  recent 
data  suggests  that  this  exponent  is  larger  than  three'"',  and  perhaps  as 
large  as  five7 .  With  this  restriction,  however,  this  class  equation  is 
useful  in  giving  a  qualitative  survey  of  some  subcritical  transition 
properties  of  liquids  and  gases. 


Table  II.  Critical  point  values  for  designated  quantities 


State  equation  numbers  correspond  to  those 
of  Table  I,  Value  of  C$  was  assumed 


3R/2 

to  compute  last 

column. 

Equation 

No. 

d2PA 

a# 

(CV  A  “CV  )/® 

(CvB-Cv)/R 

aA/(RTc)2 

2 

4 

0 

4 

0 

.612 

3 

4 

0 

4 

0 

1.633  Zc 

4 

32  f 

-6 

20  i 

4 

.563 

5 

32  f 

-6 

24 

4 

1.501  Zc 

6 

0 

8 

0 

1.622  zc 

7 

24 

0 

-1 

0 

1.871  zc 

r 


Figure  k.  Vaporizing  signal  speed  of  the  saturated  liquid  for 
equations  numbered  in  Table  I.  The  curve  shown  for 
No.  6  was  computed  for  Zc  =1/3.-  The  value  3R/2 
was  used  for  C$ . 
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APPENDIX 


I.  Showing  B  must  be  Independent  of  y 

If  B  were  allowed  to  be  a  function  of  y,  there  would  be  additional 
terms  in  Cvb(V,T)  of  Eq.  (33)*  These  extra  terms  would  be 

-  {(9/8)  Ry/ CH-B )} { 2B /  +  y(b')2/(t1-b)}. 


Since  T]  can  approach  B  in  the  single  phase  of  this  fluid,  the  last  term 
would  contribute  a  dominant  negative  value  to  CvB(V,T),  unless  B'^y) 
were  identically  zero.  Because  of  the  stability  condition  Cv(V,T)  >  0, 
it  is  then  necessary  that  B(y)  be  a  constant.  Restrictions  are  also 
required  on  A(y)  and  C(y),  but  they  are  not  as  simple.  With 
x  =  rn+CtY)]-1,  and  the  definitions  L(y)  =  2A"  (y  ) ,  M(y)  =  Zk'c'+AC" 
and  N(y)  =  2AC",  then  Eq.  (33)  becomes 


Cvb/R  =  C$/R  +  (9/8)  yx  Cl(y)  -  xM(y)  +  x?N(y)},  (A-l) 

where,  for  y  ^  1,  the  permissible  x  range  over  a  positive  interval  from 
0  to  xi  =  +  [B  +  C(y)]-1  <  +  ooj  for  y  £  1,  permissible  values  of  x  range 
over  a  similar  interval,  excluding,  however,  a  sub-interval  corres¬ 
ponding  to  7)  values  within  the  coexistence  region.  (For  this  sub-inter¬ 
val,  the  proper  expression  for  Cv  is  CVA,  given  by  Eq.  (36)}.  Eq.  (A-l) 
is  a  cubic  polynomial  in  x,  and  can  have  a  variety  of  shapes,  depending 
on  the  coefficients  L(y),  M(y)  and  N(y).  The  functions  A(y)  and  B(y) 
must  be  restricted  so  Cv  is  always  positive.  This  requirement  is 
satisfied  for  all  the  example  equations  shown  in  Table  I. 

II.  Reduced  Variables  for  the  [V-M]  Fluid 

The  set  of  values  given  in  Table  III  for  the  [V-M]  system  can  be 
used  with  the  transformation  equations  to  compare  easily  the  sub-critical 
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Table  III.  Reduced  values  of  coexistence  temperature, 
saturation  densities,  and  vapor 
pressure  for  the  Van  der  Waals- 
Maxwell  state  equation 


u 

5i 

wA 

1.000 

1.0000 

1.0000 

1.0000 

.998 

1.0902 

.91140 

.99202 

•  994 

1.1571 

.84773 

.97617 

.990 

1.2035 

.80454 

. 96048 

.98 

1.2894 

.72669 

.92191 

•  97 

1.3558 

.66844 

.88429 

•  96 

1.4121 

.62042 

.84762 

•  95 

1.4617 

.57901 

.81188 

.93 

1.5482 

. 50931 

.74318 

.90 

1.6573 

.42574 

.64700 

.85 

1.8071 

.31973 

.50449 

.80 

1.9327 

. 23967 

.38336 

.75 

2.0424 

. 17721 

. 28246 

.70 

2 . i4o4 

. 12802 

. 20046 

.65 

2.2296 

.89475xlO_1 

.13584 

.60 

2.3116 

.59778x10_1 

.86869xl0_1 

.55 

2.3875 

.37580X10-1 

.51580x10-1 

.50 

2.4585 

. 2l747xl0_1 

. 27789x10_1 

.45 

2.5251 

. 11217xlO_1 

.12134xlO_1 

.40 

2.5879 

.49109xl0“2 

.51745x10"2 

.35 

2.6475 

. I6875xl0-2 

.I5673xl0"2 

.30 

2.7042 

. 399°7xl0-3 

,31882x10"3 

.25 

2.7583 

.51259xl0“4 

.34i65x10‘4 

.20 

2.8012 

. 22296x10-6 

. Il891xl0-6 

.15 

2.8602 

.41565x10"7 

.43602x10"® 

.10 

2.9083 

.21612x10-12 

.57631X10-13 

.04 

2.9640 

.56824x10-34 

.606l2xl0'35 

0.00 

3.0000 

0.00 

0.00 
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thermodynamic  predictions  for  any  state  equation  belonging  to  this 
class.  Five  significant  figures  are  given  for  every  function  except  u, 
which  is  assumed  to  be  exact.  These  values  were  computed  at  this 
laboratory  with  the  assistance  of  D.  C.  Mylin  and  F.  H.  Macintosh.  R. 
Barieau1  gives  a  much  more  extensive  tabulation  of  various  thermodynamic 
properties  for  the  [V-M]  fluid  for  .25  ^  u  ^  1.0.  For  these  u,  the 
above  values  agree  with  his  tabulations. 
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